(N 



O 

O HOMODYNE CHARACTERIZATION 

S \ OF ACTIVE OPTICAL MEDIA 

^^ ■ G. Mauro D'Ariano, Matteo G. A. Paris and Massimiliano F. Sacchi 

c7^ , Theoretical Quantum Optics Group, INFM and Dipartimento di Fisica "A. Volta" 

Universitd di Pavia, via Bassi 6, 1-27100 Pavia, ITALY 

Abstract An effective maximura likelihood method is suggested to characterize 
the absorption/amplification properties of active optical media through 
homodyne detection. 
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^^ I The quantum characterization of optical media is an important issue 

in modern optical technology, since the noise in optical communications 
^5 , and measurements is ultimately of quantum origin. For negligible satu- 

^D I ration effects, the propagation of an optical signal in an active media is 

governed by the master equation 



Dh: 



c^ 



Q = GiL[a]g + G2L[a^]g, (1.1) 



P3 

j^ ■ where g is the density matrix describing the quantum state of the signal 

^ . mode a and L[0] denotes the Lindblad superoperator L[0]^ = O AO^ — 

5^1 ^0^0 A — i^AO'^O. If we model the propagation as the interaction of 

a traveling wave single-mode a with a system of N identical two-level 
k^ \ atoms, then the absorption Gi = jNi and amplification G2 = JN2 

V^ ' parameters are related to the number A'^i and A'^2 of atoms in the lower 

and upper level respectively. The quantity 7 is a rate of the order of 
the atomic linewidth [1], and the propagation gain (or deamplification) 
is given by G = exp[(G2 - Gi)t]. 

An active medium described by the master equation (1.1) represents 
a kind of phase-insensitive optical device. In this paper, we want to 
evaluate the parameters Gi and G2 by the maximum-likelihood (ML) 
estimation applied to data coming from random phase homodyne de- 
tection on the signal exiting the medium. The present investigation is 
motivated by the fact that ML approach has been already successfully 
applied to estimation of the whole quantum state [2] as well as to deter- 
mination of some parameters of interest in quantum optics [3] . 

Let us start by reviewing the ML approach. Let p{x\X) the probability 
density of a random variable x, conditioned to the value of the parameter 
A. The analytical form oip is known, but the true value of the parameter 



A is unknown, and should be estimated from the result of a measurement 
of X. In our case A is the couple of parameters Gi and G2, and p is the 
probability density of (random phase) homodyne data. Let xi,X2, ■■■, xn 
be a random sample of size N. The joint probability density of the 
independent random variable xi,X2, ■■■, xn (the global probability of the 
sample) is given by 

CiX)=U^^,pixk\X), (1.2) 

and is called the likelihood function of the given data sample. The maxi- 
mum-likelihood estimator of the parameter A is defined as the quantity 
Xmi that maximizes C{X) for variations of A. Since the likelihood is 
positive this is equivalent to maximize 

N 

L(A)=log/:(A) = ^logp(xfc|A) (1.3) 

k=l 

which is the so-called log-likelihood function. 

Using the Wigner representation of Eq. (1.1) one can easily solve the 
corresponding Fokker-Plank equation for the Wigner function VF(a, a*; t). 
One obtains [4] 

l^(a,«*;t)=|0exp(^- '"~f^'' W(/3,r;O) (1-4) 

with g = e^«*, 2Q = (Gi - G2), and 6^ = (Gi + G2)(l - <7^)/(4Q). The 
theoretical homodyne probability at phase (j) is simply obtained as the 
following marginal distribution 

p{x; (/))=[ dlma W{ae'^, a*e-'^) x = Re(a) . (1.5) 

For input coherent state with amplitude ao, one has W{P,(3*;0) = 
- exp(— 2|/3 — oqI^) and the convolution in Eq. (1.4) gives 

The corresponding theoretical homodyne distribution is then given by 



For non-unit quantum efficiency f] < 1, one has the replacement 

,^4_,nf + i^. (1.8) 



Homodyne characterizationof active optical media 



We applied the ML approach to determine Gi and G2 starting from 
random phase homodyne detection [(p in Eq. (1.7) randomly distributed 
in (0, tt)]. As a input reference signal we used coherent state of fixed 
known amplitude. Notice that the use of coherent states is not simply 
a matter of computational and experimental convenience. In fact, there 
is no advantage in using e.g. squeezed states, because of the phase- 
insensitive character of the device. Compare, on the contrary, the case 
of phase estimation in Ref. [3] . Some results from Monte Carlo simulated 
experiments for both the absorption (Gi > G2) and the amplification 
(Gi < G2) regime are shown in Table 1. Notice also that the case 
Gi = G2 corresponds to the estimation of Gaussian noise, since one has 
the solution of Eq. (1.1) in the form 



J vrGii 
where D{(3) = e^" -/3*a denotes the displacement operator. 



(1.9) 



Gi 


G2 


(G'i)ml 


(<5Gi)ml 


{G2)mL 


{5G2)ml 


3. 


1. 


2.97250576 


0.03489146 


0.96966708 


0.03299910 


3. 


2. 


2.93669546 


0.04629955 


1.94330199 


0.04412476 


3. 


3. 


3.03023643 


0.07376747 


3.03199992 


0.07122468 


3. 


4. 


2.98543015 


0.09926873 


3.98150430 


0.09763157 


3. 


5. 


3.16888784 


0.06556872 


5.15783291 


0.06240839 



Table 1 Maximum likelihood estimation of loss (Gi) and gain (G2) parameters of 
master equation (1.1). Input coherent state with amplitude qq = 4 has been used, 
along with A'^ = 10'' homodyne data with quantum efficiency r) — 0.6, collected after 
an effective time t = 1. SGi and 5G2 represent the statistical error of the estimation. 



In Fig. 1 we show the behavior of the statistical errors on the maxlik 
determination of the parameters as a function of the number of homo- 
dyne data and the quantum efficiency of photodetectors. The robustness 
of the method to low quantum efficiency ry is a feature of the maximum- 
likelihood technique [2, 3]. In the present case, however, it is not surpris- 
ing [see Fig. (1)], because quantum efficiency itself can be described by 
master equation (1.1) [4]. Notice the inverse square root behaviour of the 
statistical errors versus the number N of data in the sample, according 
to the central limit theorem. 
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Figure 1 Behaviour of the statistical error {5Gi)ml (circle) and [5G2)ml (triangle) 
versus number of homodyne data A^ (left, notice the scaling 5G oc jv^'^'^"^) and quan- 
tum efficiency of homodyne detector (right) in the maximum likelihood estimation of 
Gi and G2 in the master equation (1.1). Parameters (left): d = 3, G2 = 5, 77 = 0.6, 
ct{) — h, t = 1. Parameters (right): number of data N = 5 x 10'"*, Gi — 2, G2 ~ 1, 
ao = 8, t = 1. 

In conclusion, we applied the maximum-likelihood estimation approach 
to the characterization of linear active optical media through homodyne 
detection. The resulting method is efficient and provides a precise deter- 
mination of the absorption and amplification parameters of the master 
equation using small homodyne data sample. 
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